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We analyzed the localized charge dynamics in the system of N interacting single-level quantum 
dots (QDs) coupled to the continuous spectrum states in the presence of Coulomb interaction be- 
tween electrons within the dots. Different dots geometry and initial charge configurations were 
considered. The analysis was performed by means of Heisenberg equations for localized electrons 
pair correlators. 

We revealed that charge trapping takes place for a wide range of system parameters and we 
suggested the QDs geometry for experimental observations of this phenomenon. We demonstrated 
significant suppression of Coulomb correlations with the increasing of QDs number. We found the 
appearance of several time scales with the strongly different relaxation rates for a wide range of the 
Coulomb interaction values. 
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I. INTRODUCTION 



Coupled quantum dots (QDs) are recently under 
numerous experimental HjHJjl and theoretical in- 
vestigations [!|>@>@>0)H>H)llo| due to their poten- 
tial application in modern nanoscale devices dealing 
with quantum kinetics of individual localized states 

The kinetic properties 
of coupled QDs (artificial molecules) El are governed 
by the Coulomb interaction between the localized elec- 
trons and depend strongly on the dots topology, 
which determines energy levels spacing and the coupling 
rates E3 j EH ! EH • During the last decade experimen- 
tal technique gives possibility to create vertically aligned 
strongly interacting QDs with only one of them coupled 
to the continuous spectrum states 20], 21]. This so-called 
side-coupled geometry gives an opportunity to analyze 
non-stationary effects in formation of various charge and 
spin configurations in the small size structures (7|.[l0|. 
Lateral QDs are extremely tunable by means of indi- 
vidual electrical gates [22], [23]. This advantage reveals 
in the possibility of single electron localization in the 
system of several coupled dots [24| and charge states 
manipulations in the artificial molecules. Therefore lat- 
eral QDs are considered to be an ideal candidates for 
creation of an efficient charge traps. Previous studies 
demonstrated lon g-liv ed charge occupation trap states in 
a sin gle QDs [25| |26j].(27T] and single electron spin trap- 
ping [28|. Single electron trapping in the double dot 



system was performed in [29(. The temperature of the 
trapped electron was measured and tunnel coupling en- 
ergy was extracted by charge sensing measurements. A 
full configuration-interaction study on a square QD con- 
taining several electrons in the presence of an attractive 
impurity was performed in 30]. Authors demonstrated 
that the impurity changes significantly the charge den- 
sities of the two-electron QD excited states. The effect 
of correlations was revealed in the enhancement of the 
charge densities localization within the dot. QDs were 
investigated theoretically by various methods such as 
Keldysh non-equilibrium Green- function formalism [3l| . 
re-normalization group theory [321 ] , specific approach sug- 
gested by Coleman |33j |. spin-density- functional theory 
34 1 or quantum Monte- Carlo calculations pjjj . 

In this paper we consider charge relaxation in the 
system of N interacting QDs with on-site Coulomb re- 
pulsion coupled to the reservoir (continuous spectrum 
states). The analysis was performed by means of Heisen- 
berg equations for the localized electrons pair correlators. 
We demonstrated the presence of strong charge trapping 
effects for the lateral QDs geometry. We found that on- 
site Coulomb repulsion results in the significant chang- 
ing of the localized charge relaxation and leads to the 
formation of several time ranges with strongly different 
values of the relaxation rates. We also pointed out the 
significant suppression of Coulomb correlations influence 
on the localized charge relaxation with the increasing of 
dots number. 



II. THEORETICAL MODEL 
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Let us consider relaxation processes in the system of 
N identical lateral QDs which are situated in the differ- 
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ent space points and are connected only with the single- 
central QD by means of electron tunneling processes 
with the same tunneling transfer amplitudes T. We as- 
sume that the single particle level spacing in the dots is 
large than all other energy scales, so that only one spin- 
degenerate level within the QD spectrum is accessible [eq 
in the N identical dots and e in the central one). QD with 
energy level e is also connected with the continuous spec- 
trum states. Moreover we take into account Coulomb in- 
teraction between the localized electrons within the dots 
(U-in the central one and Uo in the surrounding dots). 
Hamiltonian of the system under investigation has the 
form: 



N 



H = ea+a a + ^ e„&t b ja + ^ (Tb+ a a a + Ta+b ja ) 
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where A = e — s$ is the detuning between the energy 
levels in the dots. 

System of equations ([3]) can be re- written in the com- 
pact matrix form (symbol [] means commutation): 



i^- f G= [G, A] — i(BG + GB) 
where G is the pair correlators matrix: 



(4) 



Xfc - tunneling amplitude between the single-central dot 
and continuous spectrum states. We assume T and T k to 
be independent of momentum and spin. By considering a 
constant density of states in the reservoir v§ (which is not 
a function of energy) , the tunneling coupling strength 7 
is defined as 7 = tti/qT^. 

<z+ /cLcript^/bjcr)- electrons creation/annihilation oper- 
ators in the central dot (in the N surrounding QDs). 
c^T/cfe- electrons creation/annihilation operators in the 



continuous spectrum states (k) and n a 
filling numbers in the dots. 



'{rijcr) are electron 



We'll at first analyze filling numbers relaxation pro- 
cesses in the case when on-site Coulomb repulsion is ab- 
sent in the whole system (U = Uq = 0). Let us assume 
that at the initial moment all charge density in the sys- 
tem is localized only in one of the N QDs and has the 
value ni(0) = riQ. The filling numbers time evolution 
can be analyzed by means of kinetic equations for bi- 
linear combinations of Heisenberg operators a+ ja a and 



(t);a+a a = G aa (t) 



b Pj'« = G jf 
b+ a a<j = G ja (t);a+b jrT = G a j{t) 



(2) 



Localized charge time evolution can be obtained from 
the system of equations for the Green functions G y , 
G a j, Gj a and G aa : 



G = 



( Gaa G a l 
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and matrix A has the following form 
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The tunneling coupling matrix B has only one nonzero 
element ||B||n = 7. 

The formal solution of the system ((4]) can be found 
with the help of evolution operator: 



G(t) = e 



— J-Bt+iAt] 



G(0)e 



[Bt-iAt] 



(7) 



Consequently the average value of filling numbers time 
evolution in the one of the N QDs can be found from the 
following expression: 



( nj (t)) = Gjj{t) = 5>- St +^] JQ G ab (0)[ f 



Bt-iAt- 



hi (8) 
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Due to the condition that initial charge is localized only 
in one QD with number j, the following initial conditions 
are fulfilled: < n Jcr (0) >= Gjj(0) = n , < n aa (0) >= 
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0, < n /CT (0) >= 0, if j ^ j' and G - (0) = G aj (0) = 
G ja (0) = 0. 

Let us analyze filling numbers time evolution in the 
central QD and in the dot with initial charge. Concerning 
initial conditions one can easily find the expressions for 
filling numbers relaxation: 



n ja (t) 



iHti 



jn 0j [e 
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where operators H = A + iB and H = A — iB are 
included. Further analysis deals with the calculation of 
matrix exponent's elements. One can easily perform this 
procedure with the help of recurrent ratio similar to the 
procedure suggested by Cummings [36j . The following 
ratios for operator H elements are fulfilled: 

( Hn )jj = H ja (H n ~ v ) a: j 

(H n ) aj = ^^ o/ (^- 1 ) /i + J ff aa (#"- 1 ) Qj 



(10) 



System of equations (ITU)) enables to get recurrent ratio 
for matrix elements (H n ) a j-. 

(H n ) aj = (A + *7)(ff- 1 )a J + N\Tf(H n - 2 ) a3 (11) 

analogous equations can be obtained for matrix ele- 

~ n 

ments (H ) a j. Consequently after some calculations one 
can get: 

J l „n | iri\ 
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Where coefficients D,D,a,a,b and b are determined as: 



(13) 
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Expanding exponents in the expression ^ in a power 

H and H series one can easily obtain the following ex- 
pressions: 
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After substituting ([14]) to equations © one gets ex- 
pressions which describe filling numbers time evolution 
in the central QD n a (t) and in the QD with the initial 
charge n,j(t) in the case when Coulomb correlations are 
neglected. 
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It is clearly evident that with the increasing of QDs 
number N initial charge noj is quite fully confined in 
the initial QD even in the presence of dissipation in the 
system due to the interaction with the reservoir. 



1 „ 

limt^^njait) = n j(l - — f 
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Simultaneously for the large number of QDs N, elec- 
tron filling numbers in the QDs reveal oscillations fre- 
quency increasing as TyN. If initial charge is localized in 
the central QD, which is coupled to the continuous spec- 
trum states, one should solve system ([3]) with the initial 
conditions: < rij a {Q) >= 0, < n aa (Q) >= G aa (0) = no, 
< n j -' o .(0) >= 0. Consequently one can get the following 
expressions for the charge time evolution: 



n a a(t) = [e tHt ]a a nOa[e lHt ]aa 

n ja (t) = [e^janoale^Uj 



(17) 



The function ip — e^ Ht ^ aa can be obtained from equa- 
tion: 



^ = i(A + ij) ■ ^ + NT ■ [e iS % a (18) 



Finally, solution will have the form: 



^yji2 e iat _ e i(A+if)t & ibt _ e »(A+*7)t 

^ = ~2D~ <y a- (A+Ty) b - (A + i 7 ) ^ ^ 19) 

where coefficients a, b and D are determined by the ex- 
pressions (|13[) . Consequently, the charge trapping effect 
is absent in this situation. 



4 



We now consider the situation when Coulomb interac- 
tion between localized electrons exists within the QDs. In 
this case it is necessary to take into account the following 
interaction part of the system Hamiltonian ([T]): 

Hint = U( Q )n aa n a -,j (20) 

where index a = a(j) and Coulomb interaction val- 
ues /7(o) correspond to the central dot (surrounding dots). 
We'll take into account Coulomb interaction by means of 
self-consistent mean field approximation [10J. It means 
that the initial energy level value e have to be substituted 
by the value e = e + U- < n aa {t) > in the final expres- 
sions for the filling numbers n acj time evolution (|15[) . So 
one should solve self-consistent system of equations. 

In the presence of Coulomb interaction system of equa- 
tions for pair correlators can be written in the compact 
matrix form: 



i— G = [G,A + d] + i(TG + GT) (21) 

where matrixes A, G are determined by expressions §5§ 
and (J6j) correspondingly, and matrixes T and C can be 
written as ||r||ij = <5.a<5ji7 and ||C||-y = SijUoGjj. 

The formal solution of the system for pair correlators 
(|2T1) can be again found with the help of evolution oper- 
ator: 

G(t) = Te [l ti {A{t ' )+5{t ' ))dt ' ] ■ G(0)T ■ e [ ~ l ^ (Mt )+6(t'))dt 



As initial charge is localized in the QD with number 
j, the initial conditions are: n^(0) = rijoSij, n a (0) = . 
Then one can obtain the expressions: 

k,k' 

G aa {t) = n 0j \n aj \ 2 (23) 

where evolution operator Q = Te^ti&it^+Cit'Vdt] 
is considered. So, one can get equations for the matrix 
elements of the evolution operator 17 



neglected, because the filling numbers amplitude is pro- 
portional to 1/N 2 . So, taking into account only Coulomb 
correlations within the dot with the initial charge, one 
can simplify the system of equations (pM)) in the follow- 
ing way: 

n aj = -(N- 1)T 2 • tt aj + i(A - 7) • (l aj + iT ■ (ijj 
(ijj = IT ■ Q a j + illonoj ■ Qjj (25) 

System of equations (f25j) can be easily solved numeri- 
cally and consequently one can analyze localized charge 
relaxation processes. 

III. CALCULATION RESULTS AND 
DISCUSSION 
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FIG. 1: Filling numbers time evolution a), in the QD with 
the initial charge and b). in the central QD for the different 
number of QDs. Parameter T/7 = 0.6 is the same for all the 
figures. 
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- iT ■ tl aj + iU n 0j ■ Ify/fy, (24) 

with initial conditions for the functions Q a j(0) = 0, 
£ljj(0) — 1, Q .(0) = 0. If we are interested in the col- 
lective effects connected with the presence of large num- 
ber of coupled QDs N, the Coulomb interaction between 
localized electrons within the initially empty dots can be 



Filling numbers time evolution within the dot with ini- 
tial charge rij(t) and within the central QD n a (t) in the 
absence of on-site Coulomb repulsion is presented on the 
FigJT] The non-resonant tunneling between the dots is 
considered ((e — £q)/7 = — 1). 

It is clearly evident that filling numbers relaxation 
changes significantly with the increasing of QDs number 
N. When initial charge is localized in one of the N QDs it 
remains confined in the initial dot even in the presence of 
relaxation processes from the central dot to the reservoir 
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for the large number of dots. When one considers two 
surrounding dots only twenty percent of charge continue 
being localized in the initial QD (FigHk). But for ten in- 
teracting QDs more then eighty percent of charge is con- 
fined in the initial dot (FigQJi). This effect can be called 
"charge trapping" and the proposed system of coupled 
QDs can be considered as a "charge trap". QDs num- 
ber N increasing also leads to the decreasing of charge 
amplitude in the central QD n a (t) for a fixed value of ra- 
tio T/7 due to the effective growth of tunneling coupling 
(FigHt). 

Typical calculation results, in the case when on-site 
Coulomb repulsion is considered only in the central QD 
where localized charge is absent at the initial time mo- 
ment, are demonstrated on the Fig[2] "Charge trapping" 
effect is clearly evident with the increasing of QDs num- 
ber even in the presence of Coulomb interaction between 
localized electrons (FigJ^J ■ 

For two QDs interacting with the central one Coulomb 
correlations strongly influence on the filling numbers re- 
laxation (Fig|3k)- A critical value of on-site Coulomb re- 
pulsion exists for a given set of system parameters which 
corresponds to the full compensation of the initial neg- 
ative detuning [10]. For the smaller values of Coulomb 
interaction, correlations lead to the increasing of relax- 
ation rate in the QD with the initial charge (Fig|3k grey 
line) in comparison with the case when Coulomb interac- 
tion is absent (Fig|2k black line), due to the decreasing 
of the initial detuning value. For the values of on-site 
Coulomb repulsion larger than the critical one, positive 
detuning occurs and filling numbers relaxation rate de- 
creases as a result of positive detuning value increasing 
(Fig J5k. black-dashed line). 

With the increasing of QDs number all the effects men- 
tioned above are still valid, but they are less pronounced 
(Fig|2b). So the role of Coulomb correlations is sup- 
pressed for the large number of QDs due to the decreasing 
of electrons occupation in the central dot. 

Let us now analyze the situation when Coulomb inter- 
action between localized electrons is taken into account 
within all the N QDs which interact with the central 
one. Calculation results are presented on the Fig(3k. and 
demonstrate " charge trapping" effect with the increasing 
of the QDs number. In the case of two QDs interacting 
with the central one Coulomb correlations reveal signif- 
icantly stronger influence on the filling numbers relax- 
ation processes in comparison with the geometry when 
five dots are considered (Fig|3b). 

Again two typical relaxation regimes were revealed. 
The first one corresponds to the decreasing of initial nega- 
tive detuning value. In this regime Coulomb correlations 
lead to the increasing of relaxation rate in the QD with 
initial charge in comparison with the case when Coulomb 
interaction is absent (Fig|3~k grey and black lines corre- 



spondingly). The second one deals with the Coulomb 
energy values large enough to compensate negative de- 
tuning and to form the positive one. In this regime filling 
numbers relaxation rate decreases as a result of positive 
detuning value increasing caused by the Coulomb inter- 
action (FigJSk grey and black-dashed lines correspond- 
ingly). 

QDs number increasing also results in the increasing 
of filling numbers oscillations frequency. Filling num- 
bers oscillations frequency for the small Coulomb val- 
ues decreases corresponding to the detuning decreasing 
and increases as a result of positive detuning formation 
(Fig 12k, c). We found the growth of charge amplitude in 
the central QD with the increasing of the dots number 
when the negative detuning value decreases and ampli- 
tude decreasing when positive detuning value increases 
(FigEb). 

When initial charge is localized in the central QD 
which is connected not only to the surrounding QDs but 
also to the continuous spectrum states " charge trapping" 
effect doesn't exist at all (FigU]). 

We now introduce the possible QDs geometry which al- 
lows to perform an experimental observations of charge 
trapping effects within the single dot. The most simple 
configuration is: N similar lateral QDs interacting only 
with the single vertically aligned dot. Single vertical dot 
is also connected to the continuous spectrum states. But 
this geometry reveals a problem of initial charge local- 
ization. Initial charge can be localized in the single QD 
in the most simple way by means of the gate voltage. So 
it is convenient to have a system with N — 1 lateral dots 
and single vertical dot with the localized charge. These 
N dots interact only with the single central vertically 
aligned QD also connected to the continuous spectrum 
states. 

To conclude, we have analyzed time evolution of the 
electron filling numbers in the system of N interacting 
QDs both in the absence and in the presence of Coulomb 
interaction between localized electrons within the dots. 
It was found that Coulomb interaction modifies the re- 
laxation rates and the character of localized charge time 
evolution. It was shown that several time ranges with 
considerably different relaxation rates arise in the sys- 
tem of coupled QDs. We demonstrated and carefully 
analyzed the presence of strong charge trapping effects 
in the proposed systems. It was found that interacting 
dots can form an effective high quality charge trap. We 
also revealed the Coulomb correlations suppression with 
the increasing of QDs number. 

The QDs geometry which allows to perform an exper- 
imental observations of charge trapping effects was sug- 
gested. 

This work was partly supported by the RFBR grants. 
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FIG. 2: Filling numbers time evolution a),c). in the QD with initial charge and b). in the central QD for the different number 
of QDs a),b). N — 2 and c). TV = 5. Coulomb interaction is taken into account in the central QD: U/j = 0-black line, 
17/7 = 10-grey line and [7/7 = 30-black-dashed line. Parameter T/7 — 0.6 is the same for all the figures. 
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FIG. 3: Filling numbers time evolution a),c). in the QD with initial charge and b). in the central QD for the different number 
of QDs a),b). TV = 2 and c). TV = 5. Coulomb interaction is taken into account in the QD with initially localized charge: 
U = 0-black line, U/y = 2-grey line and U/j — 7-black-dashed line. Parameter T/7 = 0.6 is the same for all the figures. 
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FIG. 4: Filling numbers time evolution a),c). in the central 
QD with initial charge and b),d). in one of the N surround- 
ing QDs for a),b). N = 2 and c),d). N = 10. Coulomb 
interaction is taken into account in the central QD: U/'y = 0- 
black line, (7/7 = 4-grey line and U/y = 8-black-dashed line. 
Parameter T/7 = 0.6 is the same for all the figures. 



